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Abstract 


A first-order approximation to the range and energy straggling of ion beams 
is given as a normal distribution for which the standard deviation is estimated 
from the fluctuations in energy loss events. The standard deviation is calculated 
by assuming scattering from free electrons with a long range cutoff parameter 
that depends on the mean excitation energy of the medium. The present formalism 
is derived by extrapolating Payne's formalism to low energy by systematic energy 
scaling and to greater depths of penetration by a second-order perturbation. Lim- 
ited comparisons are made with experimental data. 


Introduction 

In space radiation transport, the energy loss 
through atomic collisions is treated as averaged pro- 
cesses over the many events which occur over very 
small dimensions of most materials and is referred to 
as the “continuous slowing down approximation” 
(ref. 1). The small percent fluctuation in energy loss is 
thought to have little meaning for ions of broad energy 
spectra and especially in comparison with the many 
nuclear events for which uncertainties are still rela- 
tively large. The exception of course is in the labora- 
tory testing of potential shielding materials with nearly 
monoenergetic ion beams in which the interpretation 
of the interaction of the ion beam with shield materials 
requires a detailed description of the interaction pro- 
cess for comparison with detector responses. In addi- 
tion to the validation of physical processes, a 
theoretical model of the role of straggling is essential 
to understanding the radiobiology of ion beams as 
required in evaluation of astronaut risks that must be 
minimized at least to within some regulated level. 

Energy and range straggling received considerable 
attention with the development of accelerated ion 
beams and the associated advancement of detector 
technology. The fluctuations of signals in detector 
responses were often a confusing factor in particle 
detection, with considerable emphasis given during 
the late 1950’s and throughout the 1960’s (refs. 2 
and 3). More recently with the development of radia- 
tion therapy beams, the issue of range and energy 
straggling take on added importance because the beam 
properties near the end of the particle trajectory 
become an essential part of session planning. Gener- 
ally such issues are studied experimentally (ref. 3), but 
a workable theory would greatly enhance the under- 
standing of the radiobiology and improve therapy pro- 
tocols. With the emergence of new biomedical 


accelerators, there is increased activity in the under- 
standing of ion beam characteristics. , 

Unlike the theory of stopping power, which has 
well-founded roots in quantum theory (ref. 4), most of 
the practical methods for energy fluctuations still rely 
on a simple modification of Rutherford’s scattering 
formula. The assumption is that the energy transfer is 
like free electron scattering with a low-energy cutoff 
determined by the atomic/molecular binding proper- 
ties (refs. 4 and 5). Even then the theory is applicable 
at best to only 85 to 90 percent of the ion range and 
only at energies above a few to several MeV/nucleon. 
In this paper, an attempt is made to find a well-defined 
extrapolation procedure to overcome these limitations 
to make them practical in ion beam models. Clearly, 
extending the theory applicability would be desirable, 
but this is beyond the scope of the present task. At the 
minimum, future experiments will allow empirical 
corrections arriving at a more accurate formalism sim- 
ilar to the analysis of experimental data using a para- 
metric stopping power formalism. 

Transport Theory 

The specification of the interior environment of a 
spacecraft and evaluation of the effects on the astro- 
naut is at the heart of the space radiation protection 
problem. The Langley Research Center has been 
developing such techniques and an in-depth presenta- 
tion is given in reference 1 although considerable 
progress has been made since that publication. The 
relevant transport equations are the linear Boltzmann 
equation derived on the basis of conservation princi- 
ples (ref. 1 ) for the flux density <J )j(x,Q,E) of type j par- 
ticles as 

q • v<$>j(x,n,E) = z 1 y k (n,a,E,E') 
x 4> t (j t.Q'.F) dil’ dE' - Oj{E) <t )j(x,Q,E) ( 1 ) 


where Oj(E) and Gjj,(QQ\E 9 E f ) are the media macro- 
scopic cross sections with Oj k (Q.,Cl\E,E f ) representing 
ail those processes by which type k particles moving in 
direction Q! with energy E' produce a type j particle in 
direction Q with energy E. Note that there may be sev- 
eral reactions which produce a particular product, and 
the appropriate cross sections for equation (1) are the 
inclusive ones. The total cross section Oj(E) with the 
medium for each particle type of energy E may be 
expanded as 

CSj(E) = Cj,JE) + Oj >e \(E) + o jr (E) (2) 

where the first term refers to collision with atomic (at) 
electrons, the second term is for elastic (el) nuclear 
scattering, and the third term describes nuclear reac- 
tions (r). The microscopic cross sections and average 
energy transfer are ordered as follows: 


a. .(E) ~ ICT 16 cm 2 

(A E at ~ 10 2 eV) 

(3) 

OjJE) ~ ICT 1 9 cm 2 

(A£ e) ~ 1 0 6 eV) 

(4) 

Oj r (E) ~ ICT 24 cm 2 

(A E r ~ I0 8 eV) 

(5) 


This ordering allows flexibility in expanding solutions 
to the Boltzmann equation as a sequence of physical 
perturbative approximations. Many atomic collisions 
(~10 6 ) clearly occur in a centimeter of ordinary mat- 
ter, whereas =10^ nuclear Coulomb elastic collisions 
occur per centimeter. In distinction, nuclear reactions 
are separated by a fraction to many centimeters of con- 
densed matter depending on energy and particle type. 
Special problems arise in the perturbation approach 
for neutrons for which al (£) « 0, and the nuclear 
elastic process appears as the first-order perturbation. 

As noted in the development of equation (1), the 
cross sections appearing in the Boltzmann equation 
are the inclusive ones so that the time-independent 
fields contain no spatial (or time) correlations. How- 
ever, space- and time-correlated events are functions 
of the Fields themselves and may be evaluated once the 
fields are known. (See refs. 6 and 7.) Such correlation 
is important to the biological injury of living tissues. 
For example, the correlated release of target fragments 
in biological systems due to ion or neutron collisions 
have high probabilities of cell injury with low proba- 
bility of repair resulting in potentially large relative 


biological effectiveness (RBE) and quality factor. (See 
ref. 8.) 

The solution of equation (1) involves hundreds of 
multidimensional integral-differential equations which 
are coupled together by thousands of cross terms and 
must be solved self-consistently subject to boundary 
conditions ultimately related to the external environ- 
ment and the geometry of the astronaut’s body and/or 
a complex vehicle. To implement a solution one must 
have the available atomic and nuclear cross-section 
data, which are a major task in code development. 

Transport Coefficients 

The transport coefficients relate to the atomic/ 
molecular and nuclear processes by which the particle 
fields are modified by the presence of a material 
medium. As such, basic atomic and nuclear theories 
provide the input to the transport code database. It is 
through the nuclear processes that the particle fields of 
different radiation types are transformed from one 
type to another. The atomic/molecular interactions are 
the principal means by which the physical insult is 
delivered to biological systems in producing the chem- 
ical precursors to biological change within the cells. 
The temporal and spatial distributions of such precur- 
sors within the cell system govern the rates of diffu- 
sive and reactive processes leading to the ultimate 
biological effects. 

Atomic/Molecular Interactions 

The first-order physical perturbation to the right- 
hand side of equation (1) is the atomic/molecular cross 
sections as noted in equation (3) for which those terms 
in equation (1) are expanded about the energy 
moments S m (E) as 

S m (E)= Z-eJ Gj(E) (6) 

where Zj is based on the electronic excitation energy 
and Oj(E) is the total atomic/molecular cross section 
for delivering Zj energy to the orbital electrons 
(including discrete and continuum levels). The first 
moment (m = I) is the usual stopping power, and the 
usual continuous slowing down approximation (csda) 
is achieved by neglecting the higher order energy 
moments. The second energy moment is related to 
energy/range straggling and provides corrections to 
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the ion slowing down spectrum (ref. 1). Equation (6) is 
misleadingly simple because specification of Ej and 
Oj(E) requires a complete knowledge of the atomic/ 
molecular wave functions. A many-body local plasma 
model has been found useful in approximating the 
atomic and molecular quantities for the positive 
energy moments (ref. 1). The current stopping power 
database is derived semiempirically as the Bethe 
reduction of equation (6) in terms of mean excitation 
energies and shell corrections (ref. 1). The usual rela- 
tivistic correction and the density effect correction of 
Stemheimer are included (ref. 9). 

The passing ions are not the primary mediators of 
biological injury but rather the secondary electrons 
generated in atomic collisions which transport the 
energy lost by the passing ion to the biological 
medium. The distribution of the electrons about the 
ion path is critical to evaluation of biological injury 
(refs. 6 and 7), critical to the evaluation of shield atten- 
uation properties (ref. 8), and fundamental to dosimet- 
ric evaluation of astronaut exposure risks. Such effects 
are likewise governed by equation (1). The next physi- 
cal perturbation term is the Coulomb scattering by the 
atomic nucleus and is represented by Rutherford scat- 
tering modified by screening of the nuclear charge by 
the orbital electrons using the Thomas-Fermi distribu- 
tion for the atomic orbitals. The total nuclear Coulomb 
cross section found by integrating over the scattering 
directions is related to the radiation length. 

Nuclear Interactions 

The extent of the nuclear interaction cross-section 
database required for the transport of cosmic rays 
spans most nuclear-reaction physics from thermal 
energies to energies above tens of GeV/nucleon, 
including a large number of projectile and target mate- 
rial combinations. The types of cross sections required 
for the transport involve total yields and secondary 
energy spectra for one-dimensional transport and dou- 
ble differential cross sections in angle and energy for 
three-dimensional transport. Fortunately, neutron and 
proton cross sections have been studied at some length 
in the past. Nuclear-reaction modeling is required, 
especially for light and heavy ion projectiles, to under- 
stand the basic physical processes, and to extrapolate 
the limited, available experimental data between pro- 
jectile energies and projectile-target combinations. 


A microscopic theory for the description of 
nuclear fragmentation is being developed through the 
study of the summation of the nucleus-nucleus, 
multiple-scattering series for inclusive reactions where 
a single reaction species is considered. This approach 
originated in a theory for high-energy alpha particle 
fragmentation (ref. 10) and has been extended to 
recast the abrasion-ablation model in microscopic 
form (ref. 1 1). The microscopic theory can be shown 
(ref. 1 1) to reduce to the optical-model formulation of 
abrasion (ref. 12) which in turn reduces to the geomet- 
ric abrasion model (ref. 13). The microscopic theory 
represents a unified approach where a single formal- 
ism generates all production cross sections required 
for heavy ion transport. Previously the production of 
heavy fragments, light ions, and nucleons were treated 
separately, often with disjoint assumptions. A unified 
approach is useful because the production spectrum of 
nucleons and light ions from abrasion correlates 
directly with the formation of prefragment nuclei and 
their excitation spectra. 

The microscopic approach proceeds by formulat- 
ing the multiple-scattering series for heavy ion reac- 
tions in terms of response functions for an arbitrary 
number of particle knockouts, appropriate for inclu- 
sive reaction theory and generalized to the case of 
heavy ion abrasion dynamics (ref. 11). The reaction 
dynamics for fragmentation processes are then unified 
by the development of a single function, the multiple- 
scattering amplitude, in terms of the momentum vec- 
tors of all secondary reaction products. The reaction 
cross sections for the various secondaries are then 
found by considering the phase space for an arbitrary 
final state where there n particles are abraded from the 
projectile, leaving a projectile prefragment. The decay 
of the prefragment nuclei into the final fragment opens 
the kinematical phase space further, and this descrip- 
tion will be required for predicting the final mass 
yields as well as the momentum distribution of ablated 
nucleons or nuclei. 

The description of the development of the scatter- 
ing amplitude in terms of abrasion response functions 
has been made by using the eikonal model. The many- 
body response functions are being developed as con- 
volutions of one-body response functions with the 
shell model and a correlated Fermi gas model. The 
corrections to the eikonal theory are then well-known 
and include large angle scattering corrections and the 
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many-body effects contained in the full nuclear propa- 
gator. Ablation can then be described by well-known 
statistical and resonance theories for nominal prefrag- 
ment excitation energies with a new phenomenon pos- 
sibly occurring for extremely large values in the 
excitation energy spectrum. Recent test of the model 
(called QMSFRG) has been very encouraging for 
future database generation. 

If one replaces the quantum mechanical abrasion 
cross sections by those for nuclei represented as par- 
tially transparent uniform spheres and a semiempirical 
correction to the surface energy to correct the prefrag- 
ment excitation energy when the prefragment is far 
from equilibrium, then one obtains the semiempirical 
fragmentation model (ref. 14). This model is a highly 
efficient fragmentation database code and can repre- 
sent available experimental data, even at relatively 
low energies when Coulomb trajectory corrections are 
made (ref. 14). It is not as fundamental a code as the 
microscopic theory because it is limited by the 
semiempirical correction and by the assumption that 
nuclei are uniform spheres. 

First-Order Solution Methods 

The lowest order approximation to the Boltzmann 
equation is given in terms of the atomic collision pro- 
cesses as (ref. 1 ) 

£2 • V<fy(x,Q t £) = £ c„ )at (£+e„) 

x fytr,Q,£+e„) - o j/di (E) fy(*,£2,£) (7) 

where £ n represents the atomic/molecular excitation 
energy levels. Equation (7) is equivalent to a one- 
dimensional transport along the ray directed by £2. For 
simplicity of notation, we use a one-dimensional equa- 
tion as 

8- §j(z,E) = £ <J„ at (£+e„) 

x §j(z,E+e n ) - Oj JF) <| )j(z,E) (8) 

where the subscripts at and j are dropped in the rest of 
this paper. The boundary condition is taken as 

<|>(0,£) = 5 (E-E 0 ) (9) 

where E ( , is the initial energy. 


The solution can be written with perturbation the- 
ory as 

<|> (0) (z,E) = exp(-az) 8 (E-E 0 ) ( 1 0) 

<t> ( 1 \z,E) = c z exp(-oz) I g n 5 (E+e„-E 0 ) (II) 
2 

<t> (2) U,E) = ^7- exp(-oz) 

g n g m 5(E+e n +e m -E 0 ) (12) 

and similarly for higher order terms, where 5() is the 
Dirac delta function, and e n « E 0 has been assumed 
so that <j and g n = o n /o are evaluated at E a . The aver- 
age energy after penetration of a distance z is given by 

(E) = E (> - (e) c z (13) 

where the average excitation energy is 

(e) = Ze„g n (14) 

and the sum over n contains discrete as well as contin- 
uum terms. The standard deviation ( s ) about the mean 
energy is similarly found to be 

s 2 = ((E - (E)) 2 ) = (e 2 ) oz (15) 

with 

(e 2 ) = Ze 2 g n (16) 

Similar results can be derived for the higher moments 
of the energy distribution, which depend on atomic/ 
molecular quantities through the g n terms. Consider- 
ing the nonlinear dependence of the transported spec- 
trum on the atomic cross sections o n , it is surprising 
that the transported spectral parameters depend lin- 
early on g n . Equations (13) to (16) apply when 
E 0 » (e)oz so that the energy variations in the cross 
sections can be ignored. The expressions are easily 
generalized to deep penetration as 

<E{z)) = - (si(E(y))]dy (17) 

and similarly for the standard deviation where the 
stopping power S(E) is the first moment of the energy 
transfer given by equation (6). The degrading particle 
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energy ( E(y )) is given by the usual range energy 
relations 

r £ dE' 

R(E) = (18) 

J 0 5(£ ) 

It is clear that R(E) is the average stopping path length 
for the ions. The corresponding spectrum is taken 
herein as 


where I is the mean excitation energy usually fitted to 
experimental data. In this approximation, the effects of 
tight binding of the inner shell electrons are neglected, 
and binding effects are included only when e m j n is 
greater than zero. In this approximation, the stopping 
power is given as 

S(E) = 2A(*Z a zV(A/ ffl ) ln r 4(mM_)£j (22 ) 


<K^£) 


exp{-r£-(£(c))] 2 /[2.v(;) 2 ]} 

J2ns(z) 


(19) 


where the standard deviation s(z) is given by equation 
(15) for low penetration and its general evaluation is 
the subject of the present paper. The usual continuous 
slowing down approximation is found as s — > 0. The 
evaluation requires knowledge of the appropriate 
atomic/molecular cross sections a /r 


where N is the number of atoms/molecules per unit 
volume. Two clear limitations of equation (22) exist. 
First the stopping power will decline at low energies 
as the collision becomes adiabatic and remains posi- 
tive definite. Second, shell corrections within the ln[] 
term corrects the positive definiteness to energies on 
the order of 0.5 A MeV, which leaves only the adia- 
batic region unrepresented. In evaluation of straggling, 
we follow Payne and approximate as 


Payne’s Quasi-Free Electron Approximation 

The means of deriving a first-order approximation 
from the excitation spectrum of atoms or molecules 
are given. In practice, the appropriate cross sections 
and excitation energies are not known, and the details 
of the calculation with approximate wave functions 
are tedious. Using a simple approximation which is 
adequate for a first-order theory is customary (refs. 1 5 
and 16). The cross section is represented by a contin- 
uum distribution for a heavy charged particle colliding 
with free electrons for which a long-range cutoff is 
introduced as follows. The Rutherford cross section 
for an atom is given by 


do 

dz 


nZ a Z 2 e 4 (A/m) 


£e 


2 


( 20 ) 


where Z a is the number of orbital electrons of the tar- 
get atom/molecule; Z, the projectile effective charge 
(atomic number at high energies); A , the projectile 
mass; m, the electron mass; £, the projectile energy; 
and £, the energy delivered to the electron in the colli- 
sion. The range of £ is given as 

/ ^ 4 m 

£ • = — <£<£ = — E (21) 

^min ^ ~ c - c max ^ ; 


S(E) = S(E } ) 



£, 3 EA 

<E<— 1 

2 2 j 


(23) 


where E\ = 5AI/m and a = 1 - [ln(4w£'iM/)] _l . Clearly 
a global formulation will have to overcome these limi- 
tations. Payne solved the Boltzmann equation for the 
second moment s 2 (z) by using equations (20) to (23) 
and obtained 


s (Z) = 2- 


inE 


2 r- 


(3a + I )ln 


r 4mE\ 
A I 

Jj 


W)jf 

S(E„) j 


X 


1 - 


(E(z)) 


-i3a+ 1 


(24) 


where E 0 is the initial ion energy and ( E(z )) is the ion 
mean energy at penetration depth j. The spectrum at 
depth e is approximated herein as equation ( 1 9). Note 
that most experimental data are expressed as the full 
width at half maximum (FWHM) that is related to j(j) 
as x(z) = 2.355 .v(c). 


Two known limitations of Payne’s formalism are 
the limitation at low energies expressed in 
equation (23) and the failure of equation (24) as z 
approaches 85 percent of the particle mean range. 
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Because energy straggling is proportional to the ion 
energy, we show the evaluation of equation (24) as a 
ratio of x{z)IE a for incident protons of energies of 5, 
10, 20, 50, and 100 MeV in liquid oxygen in figure 1. 
Note that there is little energy dependence in this ratio 
over this broad energy range and the highest energies 
nearly collapse into a single curve depending only on 
the fraction of penetration depth zJR (r Evidence of the 
failure of equation (24) is seen as the ever rapidly 
broadening of the transmitted energy spectrum past 
the 80-percent penetration point. To reach a global 
formalism we must resolve both limitations. 



Figure 1. Energy scaled full width at half maximum 
(FWHM) of energy distribution of monoenergelic proton 
beam as function of scaled range in liquid oxygen. 


where primed quantities are defined at the range z' and 
are evaluated numerically by using equation (19) with 
Payne’s result, which is equation (24). To the extent 
that the energy widths are a very small fraction of the 
energy spectrum over 90 percent of the range, then we 
may assume that ( E'(z )) and s'(z ) can be evaluated at 
E' = (E(z')} and the resulting integral in equation (25) 
may be performed to give 


iz+z\E) 


exp -f E - ( E(z))] 2 I 1 2 [s(z') 2 + * Eo UJ) j'U) 2 ] 


J2n^s{z) 2 + x E() (z f z') s'(z) 2 




(26) 


where 


l Eo 


(z,z) 


S[(E\z))\ 

S[(E(z'))] 


The energy spread is shown in figure 2 at different 
penetration depths in aluminum. We use Payne’s orig- 
inal method as equation (24) for depths less than 
85 percent of its range. If Payne’s approach is carried 
out all the way to the full range, the width generated is 
too large as shown in the upper curve for each pair. 
The widths generated by the two subintervals accord- 
ing to equation (26) are shown as the lower curves. 
The second limitation of Payne’s formalism arises at 
the lowest energies where the quasi-free approxima- 
tion is not adequate. We have of course followed 


Global Formulation of Payne’s Method 

The spectral dependence given by equation (19) 
with the variance given by equation (24) has obvious 
limitations some of which are now resolved. Because 
equation (19) is the solution for a monoenergetic beam 
of unit intensity at the boundary, we may use linear 
superposition to estimate corrections to equation (24). 
Namely, the fluence of ions at a penetration depth z 
may be represented as 


♦ (,«'.£) - f dC ”pH £, -W--'>>] 2 ' 2 -'V>) 

J Jin s' z 


Jin 

zj£ 

J2n s'(z) 


exp {-[£'-<£ , (;))1 2 /2.v ,2 (;)} 


(25) 



Figure 2. Full width at half maximum (FWHM) of energy 
distribution of monoenergetic proton beam as function of 
scaled depth in aluminum. Upper curves near end of range 
refer to Payne’s method; lower curves are modified widths. 
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Payne’s assumptions about the lowest energies in the 
prior paragraph by the extrapolation to end of the 
mean range. We extend the lower energies by noting 
that s(z)IE () is nearly an energy-independent function 
of z/R (> by which the result at a given low energy is 
approximated by extrapolating from a higher energy 
for the converged functions by energy scaling. For this 
purpose, the I MeV/nucleon energy width 5|(z) is used 
to scale to lower energies as follows: 

E S] [zR ,//?(£)] 

s(z) = ^ (27) 

where R \ is the range of the 1 MeV/nucleon ion. 


By following the previous formalism, the spectrum at 
any depth z is written as 


<Kz, E) = 


expf-CE - E c ) 2 /2s(z) 2 ] 
J2k s(z) 


(31) 


where v(z) is taken as the value defined by the previ- 
ous procedure for z < R () and equal to s(R 0 ) for larger 
values of z ■ The energy E c is defined as 


E 


C 


C E(Z )> (z < R a ) 

-E( A) ( A = z-R o >0 ) 


(32) 


There is also a conceptual issue to resolve. The 
energy spectrum at penetration to the mean range 
(z = R a ) is given by 


WR a ,E) 


ex P r-£ 2 /2.v(/? 0 ) 2 

s(R 0 ) 


(28) 


Equations (30) to (32) then result in 

I +erf \EJs(z ) ] 
= T i 


(33) 


The total fluence is shown for protons in aluminum at 
energies 5, 50, and 100 MeV in figure 3. 


where by definition only half the particles penetrate. 
Those particles of E > 0 will continue to penetrate past 
the mean range. The mean range to stopping is used to 
estimate the spectrum to larger depths of penetration 
as follows. For values of z = R (} + A, the effects on the 
spectrum can be estimated by using (£( A)) as 


4>(A +R 0 ,E) 


expj— [ + (E(A))] 2 /2s(R 0 ) 2 

J2n s(R o ) 


(29) 


which will vanish as A becomes large. The resulting 
total fluence will display an approximate error func- 
tion dependence with a nearly symmetric decline to 
zero in the neighborhood of the mean range. 



-.2 1 i i- u i j 1 i ■ 

.98 .99 1.00 KOI 1.02 

x/R () 


Figure 3. Total particle fluence of proton beam in alu- 
minum near maximum range as function of beam energy. 


Total Particle Fluence 

The total particle fluence as a function of depth is 
the means by which ion range is defined. The total 
particle fluence <I>(z) is given by 

0(z) = I <p(z,E) clE 


Comparison With Experiment 

Numerical results based on this simple procedure 
have been obtained and compared with limited 
experimental data published a long time ago. It is well 
known that all previous models predicted fairly well 
in comparison with experimental data when the 


(30) 


7 


penetration is well below 85 percent of the range. 
However, beyond 85 percent of its range, the models 
tend to predict a much larger width. This deficiency is 
corrected by using equation (26) and calculating the 
Gaussian width at the penetration well beyond 85 per- 
cent of its range. Namely, we approximate the width 
by equation (24) for z < 0.85 R a and extrapolate to 
larger z with equation (26). The calculated FWHM 
values and the published experimental values of 
Tschalar and Maccabee (ref. 3) are given in tables 1 
to 3. The experimental values are simply read off the 
graph from reference 3. With these limited compari- 
sons, tables 1 to 3 seem to establish that our modified 
model is capable of computing the energy straggling 
width well within experimental data at the depth cov- 
ering all the distances of penetration up to the maxi- 
mum range. 

Table 1. Calculated and Experimental Values of FWHM at 
Different Thickness of Penetration of Proton Beam of 
49.10 MeV Energy on Aluminum Target 

[R„ = 2.8283 g/cm 2 ] 


Thickness, g/cm" 

FWHM, MeV 

Calculated 

Experimental 
(ref. 3) 

2.605 

2.51 

2.45 

2.675 

2.74 

3.00 

2.713 

2.92 

3.30 

2.745 

3.14 

3.75 

2.760 

3.29 

4.00 

2.785 

3.64 

4.30 

2.802 

4.08 

4.55 

2.820 

5.25 

4.60 


Concluding Remarks 

Although the final approach to the straggling 
problem adopted in the present paper needs further 


Table 2. Calculated and Experimental Values of FWHM at 
Different Thickness of Penetration of Proton Beam of 
19.68 MeV Energy on Aluminum Target 

[R () = 0.55677 g/cm 2 ] 


'j 

Thickness, g/cirr 

FWHM, MeV 

Calculated 

Experimental 
(ref. 3) 

0.099 

0.22 

0.26 

0.267 

0.40 

0.45 

0.398 

0.60 

0.62 

0.497 

0.99 

1J0 


Table 3. Calculated and Experimental Values of FWHM at 
Different Thickness of Penetration of Alpha Particle of 
79.8 MeV Energy on Aluminum Target 

\R ( , = 0.57023 g/cm 2 ] 


Thickness, g/cm 2 

FWHM, MeV 

Calculated 

Experimental 
(ref. 3) 

0.099 

0.43 

0.50 

0.267 

0.80 

1.00 

0.398 

1.17 

1.10 

0.497 

1.84 

2.07 

0.553 

• 2.73 

2.70 


development, it does provide a usable formalism 
which seems to Fit the limited experimental data. 
There is no reason to believe that other ions are not at 
least as well represented by the present formalism. At 
least first-order estimates of the ion stopping problem 
is in hand, and an improved estimate of the near end of 
range effects on ion transport is now available in a 
simple computational algorithm. 
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